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We consider a quantum superposition of Bose-Einstein condensates in two immiscible 
internal states. A decoherence rate for the resulting Schrodinger cat is calculated and 
shown to be a significant threat to this macroscopic quantum superposition of BEC's. 
An experimental scenario is outlined where the decoherence rate due to the thermal 
cloud is dramatically reduced thanks to trap engineering and "symmetrization" of the 
environment which allow the Schrodinger cat to be an approximate pointer state. 
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<^ ■ I. MOTIVATION AND SUMMARY OF RESULTS 

*o ■ 

Microscopic quantum superpositions are an everyday physicist's experience. Macroscopic quantum superpositions, 
£NJ , despite nearly a century of experimentation with quantum mechanics, are still encountered only very rarely. Fast 
t> ■ decoherence of macroscopically distinct states is to be blamed 0|. In spite of that, recent years were a witness to 
an interesting quantum optics experiment Q on decoherence of a few photon superpositions. Moreover, matter-wave 
interference in fullerene Ceo has been observed ||. Another recent experiment has succeeded in "engineering" the 
environment in the context of trapped ions [Q. This success plus rapid progress in Bose-Einstein condensation (BEC) 
of alkali metal atomic vapors Q tempt one to push similar investigations even further into the macroscopic domain. 
The condensates can contain up to 10 atoms in the same quantum state. What is more, it is possible to prepare 
condensates in two different internal states of the atoms. Some of these pairs of internal levels are immiscible, and 
their condensates tend to phase separate into distinct domains with definite internal states The immiscibility 
seems to be a prerequisite to prepare a quantum superposition in which all atoms are in one or the other internal 
state, = (|iV, 0) + \0,N))/y/2, where N is the total number of condensed atoms. There are at least two theoretical 
proposals how to prepare a macroscopic quantum superposition in this framework Neither of them addresses 

the crucial question of decoherence. 

We find the quantum superposition state involving significant number of BEC atoms to be practically impossible 
£j . in the standard harmonic trap. Our master equation (derived in Section IV) , when applied to the standard harmonic 
trap of frequency uj, gives for the decoherence rate due to the environment of noncondensed atoms 
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*dec ~ 167F W-f VT N 1 , (1) 



where N is a number of condensed atoms, vt = \j2k^,T ' /m is a thermal velocity in the noncondensed thermal cloud, 
a is a scattering length, V is a volume of the trap, and Ne is a number of atoms in the thermal cloud, 

N E « . (2) 

Here /i < is a chemical potential. 

Equation (|[) gives the rate of decoherence to the leading order in the fugacity z = exp((3/j,) and also to the leading 
order in the condensate size, N. Next-to-leading order terms are 0{z 2 ) and O(N), in agreement with Refs. |) 10|, 



they were neglected here. We also assume that N is small enough for a condensate to live in a single mode as opposed 
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to the large N Thomas- Fermi limit. Given all these assumptions, in the standard harmonic trap setting Eq.(|lj) is a 
lower estimate for the decoherence rate. 

Even without going into details of our derivation, which are given in Section IV, it is easy to understand where 
a formula like Eq.(Q) comes from. TV 2 is the main factor which makes the decoherence rate large. It comes from 
the master equation of the Bloch-Lindblad form p ~ — [Na — Ab, [Na — Nb,p\], with A and B the two internal 
states of the atoms. TV 2 is the distance squared between macroscopically different components of the superposition 
(|iV, 0) + \0, N))/\/2 - the common wisdom reason why macroscopic objects are classical [Q. The factor in brackets 
is a scattering rate of a condensate atom on noncondensate atoms - the very process by which the thermal cloud 
environment learns the quantum state of the system. 

Let us estimate the decoherence rate for a set of typical parameters: T = l^K, w — 50Hz, and a — 3 . . . 5nm. 
The thermal velocity is « 10 _2 m/s. The volume of the trap can be approximated by V — 47ra^ oturn /3, where 
fli-Gturn = \/ 2k^T / mid 2 is a return point in a harmonic trap at the energy of k^T. We estimate the decoherence time 
as i dec « 10 5 sec/(A E A 2 ). For jV B = 10° . . . 10 4 and N = 10 . . . 10 7 it can range from 1000s down to 10~ 13 s. For 
N = 10 our (over-)estimate for idec is large. However, already for N — 1000 and Nr = 10 (which are still below 
the Thomas-Fermi regime) idec shrinks down to milliseconds. Given that our tdec is an upper estimate and that big 
condensates are more interesting as Schrodinger cats, it is clear that for the sake of cat's longevity, one must go 
beyond the standard harmonic trap setting. 

From Eqs.(|l],||) it is obvious that the decoherence rate depends a lot on temperature and on chemical potential. The 
two factors influence strongly both N~e and vt- Both can be improved by the following scenario, which is a combination 
of present day experimental techniques. In the experiment of Ref. [nTl a narrow optical dip was superposed at the 
bottom of a wide magnetic trap. The parameters of the dip were tuned so that it had just one bound state. The 
gap between this single condensate mode and the first excited state was 1.5/iK, which at T = 1/iK gives a fugacity of 
z = exp(— 1.5). We need the gap so that we can use the single mode approximation. At low temperatures, the gap 
results in a small fugacity, which is convenient for calculations. We propose to prepare a condensate inside a similar 
combination of a wide magnetic and a narrow optical trap (or more generally: a wide well plus a narrow dip with a 
single bound mode) and then to open the magnetic trap and let the noncondensed atoms disperse. The aim is to get 
rid of the thermal cloud as much as possible. A similar technique was used in the experiment of Ref. fl^| . 

Let us estimate the ultimate limit for the efficiency of this technique. At the typical initial temperature of 1/iK the 
thermal velocity of atoms is 10 _2 m/s. An atom with this velocity can cross a 1/zm dip in 10 _4 s. If we wait for, say, Is 
after opening the wide trap, then all atoms with velocities above 10 _6 m/s will disperse away from the dip. A thermal 
velocity of 10~ 6 m/s corresponds to the temperature of 10 _8 /L(K. As the factor NeVt ~ T 7 / 2 in Eq. ([!]), then already 
Is after opening the wide trap the decoherence rate due to non-condensed atoms is reduced by a factor of 10 -28 ! 

It is not realistic to expect such a "cosmological" reduction factor. The "dip" which is left after the wide harmonic 
trap is gone could be, for example, a superposition of an ideal dip plus a wide shallow well (which was a negligible 
perturbation in presence of the wide harmonic trap). The well would have a band of width AE of bound states which 
would not disperse but preserve their occupation numbers from before the opening of the wide trap. They would stay 
in contact with the condensate and continue to "monitor" its quantum state. Even if such a truncated environment 
happens to be already relatively harmless, there are means to do better than that. 

Further reduction of the decoherence rate can be achieved by "symmetrization" of the environmental state. Perfect 
symmetrization can be obtained provided that: 

i) Atoms have two internal states \A) and \B). 

ii) | A) and \B) feel the same trap potential. 

iii) A — A and B — B scattering lengths are the same. 

iv) A Hamiltonian has a term which drives coherent transitions A <-> B with a frequency of X/h. 

v) AE <C A. 

The term in (iv) can be realized by driven coherent transitions like in the experiments of Refs. |D| . 

Given an ideal symmetry between A and B (assumptions i— iii) , the eigenmodes of the dilute environment have 
annihilation operators S s ~ (a s + b s ) and O s ~ (a s — b s ), which are symmetric and antisymmetric respectively. a s and 
b s are annihilation operators for the two internal states of an atom in the trap eigenstate numbered as "s" with trap 
eigenenergy of e s . S s has energy (e s — A) and O s has energy (e s + A). The symmetric and antisymmetric A_E-bands 
of states can be visualized as two ladders shifted with respect to each other by 2A. In other words, the two sets 
of states feel the same, but shifted, trapping potentials. If (v) is satisfied, then the antisymmetric O s 's are nearly 
empty since they can evaporate into symmetric states and then leave the trap. The symmetric S a , s cannot distinguish 
between A and B so they do not destroy the quantum coherence between the Schrodinger cat's components. After 
symmetrization, Ae in Eq. (II]) has to be replaced by the final number of atoms in the antisymmetric states only; 




n x exp((/i - A)/fc B T), for 2A < AE 
0, for 2A > AE 



(3) 
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Here T is the temperature before opening the wide trap. n\ is the number of antisymmetric bound states which 
remain within the A_E-band of symmetric states. Atoms in these antisymmetric bound states cannot disperse away. 
For 2A > AE this number n\ is zero and there is no decoherence from the thermal cloud. 

In this perfect symmetrization limit the states |±) = (\N, 0) ± \0,N})/^/2 exist within a decoherence-free pointer 
subspace of the Hilbert space, since they have degenerate eigenvalues of the interaction hamiltonian V [jl4],|l5],[l6| . Any 
state of that subspace can be written as a\N, 0) + /3|0, N), with a and (3 complex numbers. If 'P[ q |jv.o)+/3|o,a')] denotes 
a projector onto that subspace, then 

[y,"P[a\N,o)+0\o,N)]] = 0, (4) 

which means that any quantum superposition a\N, 0) +/3|0, N) in the subspace is an eigenstate of the system operators 
in the interaction hamiltonian (a perfect pointer state), and as such will retain its phase coherence and last forever. 
The interaction Hamiltonian between the condensate and the thermal cloud is a sum of products of condensate 
operators and environmental operators. Only terms with symmetric environmental operators are relevant because 
the antisymmetric states are empty. The total Hamiltonian is symmetric with respect to A <-> B so, to preserve this 
symmetry, the relevant terms with symmetric environmental operators also contain symmetric condensate operators. 
The argument simplifies a lot for small fugacity where there is only one leading term with the Na + Nb condensate 
operator. The states |±) are its eigenstates with the same eigenvalue N. They are also (almost) degenerate eigenstates 
of the condensate Hamiltonian build out of Na,b- The coherent transitions A <-* B break this degeneracy of |±) but 
the difference of their eigenfrequencies is negligible as compared to the usual condensate lifetime of ~ 10s. In the 
next-to-leading order in fugacity there are symmetric interaction terms which change the number of condensed atoms. 
These terms drive the \N, 0) and |0,7V) states into slightly "squeezed-like" states \S, 0) and |0, S) respectively [ p"7| . 
There are also terms which exchange A with B. They give each state a small admixture of the opposite component. 
Superpositions of these are still decoherence-free pointer subspaces - there are no relevant antisymmetric operators to 
destroy their quantum coherence. When the antisymmetric environmental states begin to be occupied (see Eq.(||)), 
then the commutation relation Eq.(^j) is only approximate and states within the subspace will decohere. To leading 
order in fugacity and condensate size, the decoherence rate is given by Eqs. ([fj) and (^). 

This pointer subspace is not perfect - its existence is in apparent contradiction with the finite lifetimes of the 
condensates which can last for at most 10 . . . 20s. The reason is that the thermal cloud is not the only source of 
decoherence. The condensate loses atoms because of Rayleigh scattering, external heating, and three-body decay. 
The atoms which escape from the condensate carry information about its quantum state. They destroy its quantum 
coherence. Three body decay, the last process of the three above is the most important one fLl|| . In the experiment 
of Ref. [|ll| the measured loss rate per atom was 4/s for N = 10 7 or around 1 atom per 10 _ ^TThe last rate scales 
like N 3 so already for N = 10 4 just one atom is lost per second; decoherence time is Is. Another possibility is to 
increase slightly the dip radius. The loss rate scales like density squared so an increase in the dip width by a factor 
of 2 reduces the loss rate by a factor of 2 6 = 64. 

Ambient magnetic fields are yet another source of coherence loss [H. The condensed atoms have magnetic mo- 
ments. If the magnetic moments of A and B were different the magnetic field would distinguish between them and 
would introduce an unknown phase into the quantum superposition, thus rendering its underlying coherent nature 
undetectable. Fortunately the much used \F,mp) = |2, 1), |1, — 1) states of 87 Rb have the same magnetic moments. 
For them the magnetic field is a "symmetric" environment [jl8|. 

One more source of decoherence is the typical » 1% difference between the A — A and B — B scattering lengths 
which violates the assumption (iii) above. The Hamiltonian is not perfectly symmetric under A <-> B. Even for 
a perfectly symmetrized environment symmetric environmental operators couple to not fully symmetric condensate 
operators. This means that for the 1% difference of scattering lengths symmetrization can improve decoherence time 
by at most two orders of magnitude as compared to the unsymmetrized environment. 

In summary, we outlined a BEC scenario for an experimental realisation of a decoherence-free Schrodinger cat. 
This scenario has two ingredients: 

1) opening of the wide trap followed by an evaporation of the thermal cloud, 

2) and symmetrization of the environment. 

The Schrodinger cat is expected to be a quantum superposition of number eigenstates, \N, 0) and |0,iV). 

More details can be found in Sections II-IV, where we study the decoherence rate and the idea of symmetrized 
environment. Finally, Section V contains discussions. 

II. SYMMETRIZED ENVIRONMENT 

We introduce annihilation operators for \A) and \B), 
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4>a{x) = a g(x) + ^ a s u s (x) 

S 

(j) B {x) = b g(x) + ^ b sU s (x) 



(5) 



where g{x) is the ground state wave- function localized in the dip with energy e g < 0, and s is an index running over 
excited (environmental) states with ortonormal wave- functions u s (x) . The Hamiltonian of the system is 



H = d 



v^UaK^Ub) + {^AA^a + V0 + A V^ + U(r)^ A <f> A ~ MUb} + {A <-► B} 



(6) 



Here u = AttU ciAA/m {o-aa = o-bb are the inter-scattering lengths) and v = Anh aAs/m (o,ab is the intra-scattering 
length). The immiscibility assumption implies v > u. U{r) is the trap potential and A is the strength of the coherent 
driving transitions A «-> B. Substitution of (|^) into (|(]) and subsequent linearization gives a Hamiltonian for the dilute 
environment 



He = ^2 M a K + b t b s) - H a % + bla 3 )] 



(7) 



e s is the single particle energy of level s. We take the lowest enviromental energy min[e s ] — 0. A transformation 

a s + b x _ a s - 6 S 



S s 



V2 



V2 



brings He to a diagonal form 



(8) 



(9) 



Symmetric S"s and antisymmetric O's form two identical ladders of states but shifted with respect to each other by 
2A. There is a gap (—A — e g ) > between the lowest S s state and the ground state. In equilibrium the state S s 's 
occupation number is nf = l/[exp(/3(e s — A — e g )) — 1] w zexp(— /3(e s — A)), where the fugacity z = exp(— /3|e ff |) is 
assumed to be small. The O s 's occupation number is n® « rif exp(— 2/3A). If 2A/3 3> 1, then O's would be nearly 
empty. The condensate two-mode Hamiltonian is 



H c = e s (a f a + 6 f 6) - A(a f 6 + 6 T <z) + y (aVaa + 6+6+66) + Wc ( a tfot a6 ) > 



(10) 



where e.g. u c — u J d 3 x g 4 (x) and u c > u c . This Hamiltonian was studied in detail in 0,^9). For a purity factor 
e = (A/(u c — u c )N) <C 1 the lowest energy subspace contains two macroscopic superpositions 



|+) = ^=[(a + )* + 10,0) = ^=(\N,0) + |0,JV» 

R = j=Ko!) N (b r ) N ] 10,0) = ^=QN,0) - |0,JV» 



(11) 



The lower |+) and the higher |— ) states are separated by a small energy gap of N(u c — w c )elne. If we just have 
e < 1, the |±) states contain an admixture of intermediate states \N — 1, 1), . . . , |1, N — 1) such that their overlap is 
(+|— ) = e. For e » 1 they shrink to 



(a* ±6+)* |0,0) 



(12) 



^From now on we assume the pure case, eCl. 

Finally, the interaction Hamiltonian V, which contains all condensate-noncondensate vertices, is a sum of a sym- 
metric 



V s = [(4u + 2w)(a t a + 6 f 6) + (2«)(a t 6 + 6 f a)] 



[(u + w)(a t +6 f )] <g> 



h.c. 



4 



[(4u + 2v)(a t a + tfb) - {2v){o)b + tfa)] 
[at + &t] 



E °l°^s P 



+ 



S1S2S3 S1S2S3 



+ h.c. 



(13) 



and an antisymmetric part 

V = [(4it - 2u)(a t a - tfb) + (2u)(a f 6 - 6 f a)] 



+ 



.S1S2S3 

The coefficients are given by integrals 

1 



E Sl ^"^ + h - c 
S11S2 

E Af .,.,((« + v)OlO S2 S3 + (« - v)0lS n S n + (2u)SlS S2 S3 ) 



+ h.c. 



d 3 a; g 2 u* Sl u S2 , /3 s j s * S3 



2%/2 



d 3 a; gu* Sl u* 2 u s 



(14) 



(15) 



In Vs and Vb we neglected vertices C + C — > C + NC and C + C — > NC + NC, where C is a condensate and NC is 
a noncondensate particle, and their hcrmitian conjugates. They are forbiden by energy conservation due to the gap 
between the condensate mode and the lowest environmental state. 

All the terms in Vs,o were arranged in the form [condensate operator]® [environment operator]. Vs contains only 
[c.o.J's which are symmetric under a <-> b. They act in precisely the same way on both components of the superposition 
and as such they do not destroy coherence between the macroscopic components. To illustrate this let us calculate a 
commutator of the leading order term in Vs with the projector onto the subspace of states a\N, 0) + (3\0, N), 



(4u + 2w)(a t a + 6 t 6)(g) E Q S S2 a s * S2 , V 



a\N,0)+(3\0,N)] 







(16) 



This commutator vanishes because states of the form a\N, 0) + j3\0, N) are eigenstates of Nq = a^a + b'b. Therefore, 
this subspace would be a pointer subspace fl4|| if this leading term were the only term in the interaction Hamiltonian. 
However, Vs has other terms that are not simply functions of the total condensate number operator. Coherent states 
of an annihilation operator a + b are exact eigenstates of the [c.o.] a + b and approximate eigenstates of a' +6'. 
These coherent states, however, are combinations of ( |l2| ) and as such they are not in the lowest energy subspace of 
Hq. What is more, the matrix elements of a + b are of the order 0(\/~N) which is negligible as compared to the 
matrix elements of the number operator. If we project a + b and its h.c. on the subspace (11), then their approximate 
eigenstates for large \z\ are coherent states a\z, 0) + f3\0, z). The decoherence effects oi a + b and Nq put together lead 
to a superposition of macroscopic "squeezed-like" states a\S, 0) + /3|0, S). This result is by now well established for a 
single component condensate ]1^]. Finally, the [c.o.] a^b + b^a drives the state out of the pointer subspace (|li"|). Its 
effect is supressed by the purity condition e«l and is also negligible as compared to the direct effect of the A-term 
in He- Similar comment applies to the "out of the subspace" action of + fot . Therefore, the effect of all these terms 
in Vs imply that the subspace spanned by \N, 0) and |0,iV) is no longer an exact decoherence- free pointer subspace, 
i.e. [Vs, V[ a \N,o)+i3\o,N)]] 7^ 0- The correct pointer subspace would be one spanned by those "squeezed-like" states. 
However, to leading order in fugacity and in the condensate size, TV, the dominant terms are those depending on iVc; 
the subspace of a\N,0) + f3\0,N) is (approximately) the exact decoherence-free pointer subspace. 

The antisymmetric [c.o.J's in Vb act in opposite way on both components; they destroy their quantum coherence. To 
illustrate this let us calculate a commutator of the leading term in Vb with the projector operator onto the subspace 
a\N,O)+0\O,N), 



|JV,0)+/3|0,iV>] 



(4u - 2v)(a^a - 6+6) ® f E ^ S2 a sJ S2 + h.c. j , V [a] 

( E SlO S2 a slS2 +h.c.\ 

\S1,S2 / 



2N(Au - 2v) (a@*\N, 0)(0, N\ - a*/3\0, N)(N, 0|) 



(17) 



The commutator is O(N). The occupation numbers of O's are supressed by a Boltzmann factor n exp(— 2A/3); for 
2X(3 3> 1 the O's are unoccupied. As such they cannot scatter into condensate particles. It is also impossible to 
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scatter S's into O's thanks to energy conservation. The unoccupied O's are irrelevant for decoherence and as such 
can be neglected in V and in the above commutator. This effectively sets the dangerous Vo to zero and leaves us 
only with Vs 7^ 0. fn the absence of O's the states a\N, 0) + (3\0,N) are a pointer subspace (or a decoherence-free 
subspace). The symmetrized environment of just S"s defines the quantum states of the components but it does not 
destroy their mutual coherence. 



III. MASTER EQUATION 

Symmetrized environment is a robust idea whose validity does not depend on detailed calculations. Nevertheless, 
for the sake of illustration we derived (by a perturbative expansion in V) an approximate Bloch-Lindblad form 
master equation for the reduced density matrix p(t) of the condensate modes. The calculations are long but rather 
straightforward, their details can be found in Section IV. Here we just give the final result. 

p= l -[p,H^]+h+Po ■ (18) 
H™ n is a renormalized condensate Hamiltonian, 

Hg n = H C + (H c ) + ciAj Ai + c 2 A^A 2 + cgA^Ag + c 4 A|A 4 . (19) 

(. . .) means an average over an initial environment thermal density matrix and p$ and po are contributions from Vs 
and Vo respectively. They read 

p s = diD[Ai]p + d 2 D[A 2 ]p , (20) 
Po = d 3 D[A 3 ]p + d 3 D[Al]p + d 4 D[A 4 ]p (21) 

with D[A]p = {A^Ajp} — 2ApAt the Lindblad operator. The operators introduced in this expression are defined as 

Ai = (4it + 2w)(a t a + & t 6) + (2w)(a t 6 + 6 t a) , (22) 

A 2 = (u + «)((»* +6 + ) , (23) 

A 3 = (4u-2v)(a t a-& t &) + (2t;)(a t 6-6 t a) , (24) 

A 4 = a) - 6 f . (25) 

The coefficients c, and di are proportional to fugacity z. What is more 03,4/01,2 ~ e~ 2A/3 and d 3 ^jd\^ ~ e~ 2A/3 . 
As anticipated, the antisymmetric contributions are supressed by a Boltzmann factor due to O's. This approximate 
master equation confirms our heuristic arguments. 

Any asymmetry between A and B is a source of decoherence. If, for example, the A — A and B — B scattering 
lengths were different, uaa = u + Su and ubb — Su, then Vo would acquire extra terms <~ 5u. They would show up in 
po with coefficients d ~ z u8u. The ratio of these to the "symmetric" coefficients d\, 2 is Su/u (which is approximately 
0.03 for \F — l;rriF = —1) and \F = 2]mp = 1) states of 87 Rb ) as compared to d 3 ^/di^ ~ exp(— 2A/3). If Su/u is 
less than the Boltzmann factor, then this asymmetry is not a leading source of decoherence. 



IV. DERIVATION OF THE MASTER EQUATION 

We derive a master equation by perturbative expansion in V or in the coupling constants u, v which we regard to 
be of the same order. We assume the initial density matrix at t — to be a product p(0) — p(0) ® Pe of the system 
(condensate) and environment (noncondensate) density matrices. To take a more accurate starting point we make a 
rearrangement V — > V — (V) and He — > He + (V) where (■•■)= Tte[. • • Pe] is a trace over the environment thermal 
density matrix at the initial time. The new Hq differs from the old one by a renormalization 

e g + (Au + 2v)^n s s a s * s , 

s 

A- (2v)^nfa s . s . (26) 
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We compute pit) up to second order in the perturbation Hamiltonian V — J2i ® Ei, with A* an operator for the 
condensate and Ei one for the environment. The master equation for p{t) alone is obtained by tracing p(t) over the 
non-condensed modes. The result is [^0| 

P(t) = Up, H C ] f drg/^WMO), ^(-rj.p]] - -L £ /' ^^(^[^(O), {A,(-t), p}] , (27) 

ft zift Jo In 4 - Jo 

where 

9 f m (r) = ({ J E; i (r), J B i (0)}) , (28) 
3 J 7 m (r) = ([ J E i (r) !J E; j -(0)]) . (29) 

We have used that after the rearrangement V — > V — (V) and iJc — > #c + (^0) there are no linear terms in the 
perturbative expansion because (V) = 0. We first start to study the contribution of the symmetric terms Vs in the 
interaction Hamiltonian, and later we shall deal with the antisymmetric ones Vo- It is easy to show that the are no 
cross terms S-0 in the calculation of the master equation. 



A. Free Hamiltonian 

The next step in the calculation is to solve the dynamics of the free condensate Hamiltonian 

H c = (u c -v c )J 2 z +\ cS J x , (30) 
where we introduced angular momentum operators 

J x = ^b + tfa) , (31) 
J v = Utfa-<fb) , (32) 
J z = ^{aU-b%), (33) 

and discarded constant terms proportional to the total number of particles N = a^a + tfb f|. The Heisenberg 
equations of motion for the condensate operators a(t) and b(t) are 

i%a = efa - X cff b + (u c a^a + v c tfb)a , 

ihb = efb - X cff a + {u c b% + v c a)a)b . (34) 

These equations cannot be solved exactly for A cff ^ 0, when numbers of o's and b's are not conserved independently 
due to the coherent transfer of particles from one state to the other. In their mean-field version, these equations 
correspond to the well-known macroscopic self-trapping equation, which have the feature that below a critical "purity" 
the oscillations between states A and B are not complete; for e <C 1 the systems is self-locked in either of these states 
pl| . In this limit we can set A eff to zero, and for the immiscible case v c > u c the ground state of He corresponds to 
maximum eigenstates of Jf , i.e. pure cat states (|IT|). Also in this case the mean- field solutions are 

a(t) = a (0)e-W+"(«=+»c)/2] )6 ( t ) =6 ( ) e -^[ e f+Ar(«o+«c)/2] , (35 ) 
We believe that these solutions are qualitatively correct as long as the purity factor e <C 1. 

B. Symmetric Interaction Hamiltonian 

The symmetric Hamiltonian can be written as 

Vs = Ai£i + Ai^i + \a 2 E 2 + A 2 E 2 + h.c.l , (36) 
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where we have defined condensate operators 

Ai = (4m + 2u)(a t a + tfb) + (2u)(a t 6 + tfa) , 
Ai = (4m + 2v)(a t a + tfb) - {2v)(a t b + 6 f a) , 
A 2 = (u + u)(a t , 
A 2 = (u + w)(a f + 6 f ) , 

and the corresponding non-condensate operators 

Sl = Y a »l,*2 S k S *i ~ Y \a s *, s \ 2 StS s , 

Sl,S 2 s 

£ 1 = ^a sI , S2 Ot i0si _^K, !S |2 to s , 

si,s 2 s 

^2= /3 S *,S2,S 3 ('S'l 1 , 5 , S2'S's3 + Oj 1 S , S2 S3 ) 

E 2 = £ P'l,»,>,(Ot 1 0. a S aa +Sl l O a!l O at ) 



(37) 
(38) 
(39) 
(40) 

(41) 
(42) 
(43) 
(44) 



It is clear that the terms Ai and Ai and their corresponding environment operators have the same structure, and 
therefore will give the same qualitatively contribution to the master equation. In the following we shall keep only the 
Ai terms. 

The expectation values of multiple-points noncondensate operators are written in terms of the two-point functions 

(5t i (t 1 )5 S2 (t 2 ))=«5 SlS2 nf 1 e^^- A )(* 1 -* 2 ) e -<l tl - i2 l , 
(Oj 1 (ii)0 S2 (t 2 )) = 5 slS2 nO e i(- 1 +^)^-*=) e -<l t -^l , 
( J Sj 1 (t 1 )O aa (t a )) = 

(45) 

via Wick's theorem. Here nf'° = [ze^ £s±A ) — 1] _1 are Bose occupation numbers, and the 7 s 's are inverse finite 
lifetimes of the environmental states. We also expand these expectation values to leading order in fugacity z. The 
kernels are 

(Bi(r)Bi(O)) = (25i(0)iSi(T)>* - Y \a sh s 2 \ 2 nle^-^e- r ^<^ +0(z 2 ) , (46) 

si,s 2 

(E 2 (t)4(0) = (E 2 (0)E$(t))* = £ \(3 shS2 ^e^-^-^ [ ~ - - - 

Sl,S 2 ,S3 



2n? e -^< + < + ^ 



+ 



O -r(7° +7 S +7 ) 



+ 0(z 2 ) , 



(eI(t)E 2 (0)) = (eI(0)E 2 (t))=O(z 2 ) 



(47) 
(48) 



Given the kernels and the free evolution of condensate operators, we have to perform the time integrals in ( p7j ) . We 
calculate the master equation for late times t 3> 7 J ■ Hence the exponential decay of the propagators dominates 
their behavior, and we can apply the Markovian approximation replacing the upper limit in the time integration by 
infinity. This calculation results in two different contributions, one being a renormalization of the free Hamiltonian 



where 



Cl 



(-2 



Sl ,S2 

A Y IAi, 



|2„S 



SH C = ciAjAi + c 2 A*A 2 , 
- e S2 ) 



51 (ll+liy + (n-\e Sl -e S2 )) 



2 ' 



(49) 
(50) 



2n 



ft _1 (e ai - e S2 - e S3 + X + ef + N(u c + v c )/2) 



31 ill + 7f 2 + 7f 3 ) 2 + (^\e Sl - e S2 -e S3 +X + ef + N(u c + v c )/2)f 

o fT^esi - e S2 - e S3 + A - ef - N(u c + v c )/2) ' 

31 (7s , + 7f 2 + 7°) 2 + (^(^ - ^ - <*, + A - ef - N(u c + v c )/2)f 



(51) 
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The other contribution has the Lindblad form D[A]p = {At A, p} — 2ApA^ and reads 

p s = d^A^p + d 2 D[A 2 ]p , 

where 



d 



|2„S 



Sl (7f 1 +7f 2 ) 2 + (?i- 1 (e si -e 52 )) 2 ' 



S 2 ,S 3 I 



2/r 



7?, + 7f 2 + 7 S S , 



Sl (7 s S 1 +7f 2 +7f 3 ) 2 + (^ 1 (^ < 



Si C S2 C S g 

7 Sl + 7? 2 + 7 S S 3 



e S3 +X + ef + N(u c + v c )/2)y 



-+ 



81 uo 



(7 s 1 +7f 2 + 7s c 3 ) 2 + (^ i ( es 



e S3 +A + ef + 7VK + « c )/2)) 2 



C. Antisymmetric Interaction Hamiltonian 

The antisymmetric Hamiltonian is 



where now 



V = A 3 E 3 + A 4 E 4 + h.c. , 

A 3 = (4u - 2v)(a)a - + (2w)(a t 6 - &t a ) , 
A 4 = a t -& t , 

and the corresponding environment operators 

si ,s 2 

£4 = £ %^, S3 [(« + «)Oj 1 Oa a O a3 - (« - w)O^S- a S. 3 + 2wS s t i S S2 S3 ] 

si ,s 2 ,s 3 

The kernels are 

(E 3 (r)4(0)) = (E 3 (0)4(r)r = ]T ^^n^e^-^^ e^< + <^ + 0(z^ , 

Si,S 2 
Si,S 2 

(4(t)E 3 (0)) = (eI(0)E 3 (t))=O(z*) > 



(E 4 (r)El(0)) = (E 4 (0)El(r))* = £ |/^ >S2lS3 | 2 2(u + «) 2 n ? ie ^ e 'i-^- e '»- ^)-(% 1 +7? 2 +7 3 ) + 



2(u - D) 2 n°eT( £ 'i-^- ( '3+ 3A '- T (< +7a° 2 +7 s ° 3 ) 



( J Bi(0)S 4 (r))^(£;i(r)£; 4 (0))=O(z 2 ). 
The rcnormalization terms arc 



(5i? c = c 3 A|A 3 + c 4 A|A 4 



with 



C3 



= 12 £ 



ft - e S2 - 2A) 



- n: 



R ^ -e S2 +2A) 



51 (7f, + 7? 2 ) 2 + (JT 1 ^ - e S2 - 2A))2 Sl ( 7 o + 7 S)2 + {n -i {Csi €s2 + 2A))2 



C 4 = ^2 £ l&T, 

SI ,S 2 ,S 3 



, , 2 o ^ 1 ( e si _ e s 2 - e S3 - A + ef + N(u c + v c )/2) 

2{U + V ^ ^o +7 o +7 o 

si 's 2 1 's 3 



2„0 



2(u — v) n 



(7s ! + 7° + 7°) 2 + (TT 1 ^ - e S2 - e S3 - A + ef + iV K + „ c )/2))2 
7j _1 (e Sl - e S2 - e S3 + A - ef - AT(u c + v c )/2) 



+ 



81 (7s" + 7s S 2 + 7s S 3 ) 2 + (K-He Sl e S2 e S3 + 3A - ef - 7V( Uc + Wc )/2)) 2 



Finally, the Lindblad part is 



po = d 3 D[A 3 }p + d 3 D[Al}p + d 4 D[A 4 ]p , 



(68) 



where 

si.s 2 

^3 = -72 X! l as I< s 



T& + 7a° 



81 (7f 1 +7°) 2 + (^ 1 (e sl ~e S2 -2A)) 



2 ' 



7° + 7f 2 



S1,S2,S3 



1 (7f 1 +7? 2 ) 2 + (^ 1 (e sl ~e S2 +2A))2 ' 
S3 1 2 2( U + U ) 2 n s ° 1 



o + o + o 

ISl ' 'S2 1 'S3 



(7s , + 7° + 7°) 2 + - e S2 - e S3 ~ A + ef + N(u c + v c )/2)f 



2{u - vYn" 



7° + 1% + 7f 3 



(7.?, + 7f 2 + 7f 3 ) 2 + flT^x - e S2 - e S3 + 3A + ef + iV(u c + « c )/2))2 



D. Estimate of the Decoherence Rate 

To estimate the decoherence rate for our BEC quantum superposition, let us consider one of the terms in (|68|), namely 
p = lQv 2 d 3 [J Zl [J z ,p]], which arises form the A3 term. This term corresponds to elastic two-body collisions between 
condensate and non-condensate atoms, and induces phase decoherence. To this end we first make an estimation for 
the coefficient a s * tS2 that enters in d 3 . We assume that we have a big harmonic trap plus the dip located in its center. 
We take as the wave function for the ground state a gaussian g(x) = (mwdip/^fi.)) 3 / 2 exp(— miodipr 2 /Ti) : where o>di P 
is the frequency of the optical trapping potential in the dip. In principle, we should use thermal wave packets for the 
wave functions of the non-condensed particles p2[ . However, due to the strong localization of the condensate in the 
dip, we can approximate the non-condensed particles as plane waves. Then 

1 f giki-S e ~ik 2 -x -y / % _ \ 

a <>*> = 4 J rfV( %7f = ^ cxp r^ |fci " k2l j (72) 



On the other hand, the decay rate 7 is much smaller than w. This lets us replace the last factor in d 3 by 7r<5(?i (e Sl — 
e S2 )), where we have used the identity lim 7 ^o ^2+53 = tt8(x). Finally, replacing the sum over momenta by the 
continuum expression — > V J d 3 k — V J dfldkk 2 , and using the free dispersion relation = h 2 k 2 /2m, we find 

d 3 = -HIL f dhkfze-^e-^ 2 ^/ 2 ™ [ dfi 1 ^ 2 e-^t (1 " c ° SQ) (73) 

i6h J J 

where a is the angle between the two vectors k\ and fc 2 corresponding to the two scattered non-condensed particles. 
Since the condensate in the dip is very localized, its typical linear dimension being /dip = -\/?i/mu;di P , then /di p <C 



'thermal) where /thermal = y h 2 /2mk^T is the thermal de Broglie wavelength. This means that in the above expression 
the first (thermal) exponential factor dominates, and we can approximate d 3 by 



7rm(47r) 2 



u _ u , f dkl klze-^e-^ 2m = -^p E v T (74) 



where vt = y/2k&I '/m is the thermal velocity of non-condensed particles, and pe = Ne/V = 
V- 1 J2k ze~ l3X e~ f)h fc i/ 2m is the density of non-condensed particles in the antisymmetric states. 

Given the expression for d 3 , we can give an estimate for the decoherence rate. In the |J, M) basis, with J = N/2 
and J Z \J,M) = (l/2)(iV^ — Nb)\J,M), the master equation reads 

p(M, M') « 16v 2 d 3 (M - M'fp{M, M') = 167r 3 7 (M - M') 2 p(M, M') (75) 

For the quantum superposition state \N, 0) + |0, N), M — M' = N, so our final estimate for the decoherence rate is 
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Cc«16- 3 (4™ 2 fW)iV 2 (76) 

The only O(z) contribution to the amplitude decoherence can come from a term proportional to d^. However, 
under closer inspection turns out to be 0(z 2 ). d^ comes from a depletion/growth inelastic process such that a 
noncondensate particle in the initial state of s\ collides with a condensate particle and as a result they both end 
in noncondensate modes S2 and S3. For this to happen s± must have a sufficiently high energy to overcome a gap 
between the condensate mode and the environmental modes. Thus in this case n° = 0(z 2 ) and not O(z). 



V. DISCUSSION 



The aim of this paper was to discuss the "longevity" of Schrodinger cats in BEC's. We have shown that while in the 
standard traps decoherence rates are significant enough to prevent long-lived macroscopic superpositions of internal 
states of the condensate, the strategy of trap engineering and symmetrization of the environment will be able to deal 
with that issue. 

What remains to be considered is how one can generate such macroscopic quantum superposition, and how one can 
detect it. The issue of generation was already touched upon in Refs. [|]||. We have little to add to this. However, in 
the context of the Gordon and Savage proposal, it is fairly clear that the time needed to generate the cat state would 
have to be short compared to the decoherence time. If our estimates of Eq.(Q) are correct, symmetrization procedure 
appears necessary for the success of such schemes. 

Detection of Schrodinger cat states is perhaps a more challenging subject. In principle, states of the form (|7V, 0) + 
|0, iV))/\/2 have a character of GHZ states, and one could envision performing measurements analogous to those 
suggested in |^3| and carried out in where a 4-atom entangled state was studied. However, this sort of parity- 
check strategy, appropriate for TV < 10, is likely to fail when N is larger, or when (as would be the case for the 
"quasi-squeezed" states anticipated here [jl7|) N is not even well defined. 

A strong circumstantial evidence can be nevertheless obtained from two measurements. The first one would consist 
of a preparation of the cat state, and of a measurement of the internal states of the atoms. It is expected that in 
each instance all (to within the experimental error) would turn up to be in either A or B states. However, averaged 
over many runs, the number of either of these two alternatives would be approximately equal. Decoherence in which 
the environment also "monitors" the internal state of the atoms in the A versus B basis would not influence this 
prediction. We need to check separately whether the cat state was indeed coherent. To do this, one could evolve 
the system "backwards". However, this is not really necessary. For, as Gordon and Savage point out, when, in their 
scheme, we let the system evolve unitarily for more or less twice the time needed for the generation of the cat state, 
it will approximately return to the initial configuration. Thus, we can acquire strong evidence of the coherence of the 
cat provided that this unitary return to the initial configuration can be experimentally confirmed. 

These are admittedly rather vague ideas, which serve more as a "proof of principle" rather than as a blueprint for 
an experiment. Nevertheless, they may, we hope, encourage more concrete investigation of such issues with a specific 
experiment in mind. 
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